Density distributions for trapped one-dimensional spinor gases 
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We numerically evaluate the density distribution of a spin-1 bosonic condensate in its ground 
state within a modified Gross-Pitaevskii theory, which is obtained by the combination of the exact 
solution of the corresponding integrable model with the local density approximation. Our study 
reveals that atoms in the mp = Q state are almost completely suppressed for the anti-ferromagnetic 
interactions in both weakly and strongly interacting regimes, whereas all three components remain 
non-vanishing for ferromagnetic interactions. Specially, when the system is in the Tonks-Girardeau 
(TG) regime, obvious Fermi-like distribution emerges for each component. We also discuss the 
possible deviation of the spatial distribution from the Fermi-like distribution when the spin-spin 
interaction is strong enough. 
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I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates (BECs) of trapped alkali atomic clouds has stim- 
ulated active studies in many new regimes. From then 
on, BECs have become a popularly investigated platform 
for various effects of quantum many-body interaction in 
strongly correlated systems. When BECs are realized in 
optical traps rather than in magnetic ones, the spin de- 
grees of freedom in alkali atoms are liberated. As a con- 
sequence, such a system with internal degrees of freedom 
manifests very rich physics and many fascinating phe- 
nomena have been observed in spinor condensates, e.g., 
quantum entanglement of spins, spinor four-wave mixing, 
and spin domains, etc. H y, II H IE ■ 

On the other hand, in the research area of low- 
dimensional physics, the effect of dimensionality reduc- 
tion in bosonic s yste m has been investigated extensively 
H i, m m liljlljg, m and is being paid more and 
more attention [I^ [131 . There has been tremendous ex- 
perimental progress towards the realization of trapped 
one-dimensional (ID) cold atom systems [iM Il9l l20l |21| . 
Very recently, several groups have reported the obser- 
vation of a ID Tonks-Girardeau (TG) gas An 
array of ID quantum gas is obtained by tightly confining 
the particle motion in two directions to zero point oscil- 
lations [3 by means of two-dimensional optical lattice 
potentials. By loading the condensate in optical lattice 
and changing the trap intensities, and hence the atomic 
interaction strength, the atoms can be made to act ei- 
ther like a condensate or like a TG gas. The TG gas 
provides a textbook example where atom-atom interac- 
tion plays a critical role and mean-field theory fails to 
obtain reasonable results '25*1 . The dimensionless inter- 
action parameter 7 — mg/fi^ p governs the crossover from 
the weakly interacting condensate to the strongly inter- 



acting TG gas, where g is an effective ID interaction 
constant, m is the mass of the atom, and p is the ID 
density of the quantum gas. 

Within the mean field theory, Zhang and You re- 
cently studied the spin-1 atomic condensate in a cigar- 
shaped trap by solving an effective quasi-lD nonpolyno- 
mial Schrodinger equation p^ . As is well known, the 
mean field theory fails to work if the system enters the 
strongly interacting TG regime ft is therefore desir- 
able to investigate the ID spinor gases within a theoreti- 
cal approach valid to the strongly interacting TG regime 
as well as the weakly interacting Thomas-Fermi regime. 
In this paper, we investigate the ground state properties 
of the ID spinor Bose gases under the local density ap- 
proximation |9l ll4ll2^ combining with the exact result of 
the solvable Bose gas model |2^ HO] . The dependence 
of density profile on the magnetization and the effects of 
many-body interaction are studied in both the weakly 
interacting regime (7^1) and the TG regime (7 3> 1). 
The density profile for the ID spinor Bose gases in the 
ground states is then obtained by numerically solving 
the modified Gross-Pitaevskii Equations (GPEs). Par- 
ticularly, when the system is in the TG regime, we found 
that the effect of spin-spin interaction could be dramatic. 

The paper is organized as follows. In Sec. 11, we give 
a brief review of the ID spinor model trapped in an ex- 
ternal potential and derive the modified GPEs. In Sec. 
Ill, we introduce our numerical procedure for solving the 
coupled GPEs. Aiming at the realistic systems which 
may be accessible in experiments, in Sec IV, we present 
results of the ground state density distributions for the 
trapped spinor gases in different interacting regimes. A 
brief summary is given in Sec. V. 



II. FORMULATION OF THE MODEL AND 
METHOD 
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The s-wave scattering between two identical spin-1 
bosons is characterized by the total spin of two coUid- 
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ing bosons, or 2, and we denote the corresponding 
scattering lengths by ag and 02- When the cold atoms 
were trapped intensively in transverse direction with the 
transverse trapping frequency hujj_ greatly exceeding the 
chemical potential, the radial motion of atoms is es- 
sentially 'frozen out'. The dynamics is thus governed 
by a ID Hamiltonian with the effective ID interaction 
strength given by ^8, 31, 32] 



U0.2 = 



2r 



(1) 



-^0,2 



(l-C(ao,2/rf±)), 



where d± = ^JhlmiLiy^ and C 1.4603. For arbitrary hy- 
perfine states in the hyperfine manifold in which there are 
two spin-1 atoms, the effective interaction may therefore 
be written as 



Six,-Xj){UoVo + U2r2) 



(2) 



S {x. 



{co + C2Fi-Fj), 



where the operators Vq and V2 project the wave function 
of a pair of atoms into a state of total spin and 2. The 
coefficients cq and C2 are related to the effective ID in- 
teraction constants Uq and U2 through cq = ^0+^2 ^^^^ 
C2 = . The ground state of the system is ferromag- 



netic if C2 < 0, and anti-ferromagnetic if C2 > 0. F^, , Fy 
and Fz are three spin-1 matrices with the quantization 
axis taken along the z-axis direction. 

For a spin-1 Bose condensate trapped in an external 
potential Vext {x) which is independent of the internal 
states, the Hamiltonian can be expressed in the second 
quantized form as 



H 



2m dx^ 



J dx : : 



Co 
' 2 

-| / dx^tmFnhiFv)M^,^i^ 



(3) 



where {x) (4* ■ (x) ) is the field operator that annihi- 
lates (creates) an atom in the i-th internal state at lo- 
cation X, and i — -|-,0,— denotes the atomic hyperfine 
state \F = l,mF = +1,0,-1) respectively. Summation 
is assumed for repeated indices in the above Hamiltonian 
and the pair of colons denote the normal-order product. 
Within the mean field approach the properties of a spinor 
gas are determined by the following spin-dependent en- 
ergy functional 



£ = 



dx 



Zm dx 



dx 



(4) 



where p = J2i = and e (p) = cop/2. 



For a one-component Bose gas, it turns out that such 
a mean field approach works well in the weakly interact- 
ing regime but is not able to describe the density dis- 
tribution correctly in the TG regime. However, a mod- 
ified Gross-Pitaevskii theory, within which the interac- 
tion effect is properly taken into account by using the 
Lieb-Liniger solution, proves to be able to yield accu- 
rate ground-state density distribution even in the Tonks 
limit Ulillillillll. The modified Gross-Pitaevskii 
theory essentially is based on the local density approx- 
imation together with the analytical solution of the ho- 
mogenous system. Within the local density approxima- 
tion, the chemical potential of the inhomogeneous gas is 
determined by the local equilibrium condition, 

M = ^J'ho7n[p{x)] + Vext{x), 

with /ihom [p] being the chemical potential of the corre- 
sponding homogeneous system. For the one-component 
Bose gas, Phom[p] can be obtained from the solution of 
the Lieb-Liniger model. In our present problem, the spin 
exchange strength C2 is much smaller than the density- 
density interaction cq and the model Q for Vext{x) = 
and C2 = is the integrable three-component Bose 
gas model[2^ [S^- To make progress, we work in the 
scheme of modified Gross-Pitaevskii theory 0,113,01 and 
take phom [p] from the corresponding integrable model of 
three-component Bose gas. Effectively, this is equivalent 
to replacing the energy density e (p) in Eq. Q with the 
energy density of the exactly solvable three-component 
Bose model, which takes the values in the two limiting 
cases as following 



/ N / N r Co/o/2, 7 < 1, 

^(P) = ^P^(7) = {,2^2^2/6^^ 7»1. 



(5) 



The ground state energy density of the three-component 
Bose model has a similar form as its one-component cor- 
respondence, i.e., the Lieb-Liniger model, however, here 
the density p = p+ + po + p- is the total density of three 
components. We note that in Eq. (0J the kinetic energy 
term is associated with the inhomogeneity of the gas due 
to external confinement Vext (x) ^] and e (p) repre- 
sents essentially the energy density in the homogeneous 
system, while the last term the spin-spin interaction en- 
ergy in ID. 

In the weakly interacting regime (7 <C 1) the inter- 
action energy will not change the wave function greatly 
because it is negligibly small compared with the char- 
acteristic kinetic energy of an individual atom. In the 
TG regime (7 ^ 1), however, the interaction between 
atoms is so strong that the bosonic atoms behave much 
like spinless fermions. In both cases, the spin dependent 
term can be expressed in the explicit form 



P++ P-+ 2poP- + 2p+po - 2p+p_ 



(6) 
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We note that the last two terms in the above equation 
correspond to the processes that would change the spin 
states. An atom in the mp = ^ state scatters with an- 
other atom in the mp ~ —1 state, and consequently it 
produces two atoms in the mp = state or vice versa. 
Nevertheless, these processes conserve the magnetization 
of the system M ^ J dx {F} = J dx - $!$_] 

[2^133. In order to obtain the ground state from a global 
minimization of £ with the constraints on both N and 
A4, we introduce separately Lagrange multiplier B to 
conserve M and the chemical potential /i to conserve N. 
The ground state is then determined by a minimization 
of the free-energy functional T — S — /iTV — BM.. The 
dynamics of $i is governed by the coupled GPEs 



ind^o/dt 
ifid^^/dt 

with 



[H - B + C2 {p+ + po - P-)] 

[iJ + C2(p++p_)]$0 + 2c2$+$_$*, (7) 
[H + B + C2 {p- +PQ- p+)] + C2$g$;, 



H=-^^ + Ve.t {X)+F{p) 

im dx'^ 



and 



(8) 



(9) 



By numerically solving the above equations, we will de- 
termine the ground state density distributions for the ID 
spinor Bose gases trapped in an harmonic trap Vext {x) = 
^mixf^x^ both in the weakly interacting regime and in the 
TG regime. 



III. NUMERICAL METHOD 

To simplify the formalism we choose a slightly different 
notation. Taking advantage of the fact that all distances 
and energies in the calculation can be scaled in units 
of typical length and energy of the external harmonic 
potential, we introduce the standard length unit 



n 



(10) 



and rescale the spatial coordinate, the wave function, and 
the time variable as 
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FIG. 1: (Color online) The density profiles of the spin-1 
^'^Rb condensates in the ground state for the -|- (solid line), 
(dashed line) and - (dash-dot lines) components respectively, 
(a) m=0, Thomas-Fermi regime; (b) m=0.2, Thomas-Fermi 
regime; (c) m=0, Tonks regime; (d) m=0.2, Tonks regime. In 
this figure the lengthes are in units of 1.2/^m. 



The wave function is thus normalized to 
/ '^^Si^i \4'if' = 1- With the above changes, our 
coupled nonlinear Schrddinger equations become 



id(t>+/dT = H - B + g2ip+ + po - P-) (t>+ + g2(t>l4>*-, 



id(j)Q/dT 



H + 92 {p++p-) 



2^2 



(11) 



id(t)-/dT = H + B + g2{p- + Pa - p+) (t>- + g2(t>l(i)*+, 

where 13 = NB/afwj and H ^ - (1/2) d'^/dx^ + x^/2 + 
F (p) with 

^(^)-{goP^ 7»1. 

The pair interaction constants are also rescaled as 

go = Nca/ahuj, 
-go = N^7ry2, 
g2 — Nc2/ahuj. 

By propagating the coupled GPEs Eq. Hll|l in imagi- 
nary time, we obtain the ground state of spin-1 BECs in 
one dimension. In each propagating step, the wave func- 
tion is normalized to conserve the atomic number. We 
ensure the conservation of magnetization Ai by adjusting 
the Lagrange multiplier B. In our procedure the Crank- 
Nicholson scheme is used to discretize Eq. |[TT|) [37| . We 
take the initial wave function to be a complex Gaussian 
with a constant velocity: exp[—x^/2q — ikx]. Here q 
and k are adjustable parameters that shall not affect the 
final converged ground state (361] . 



4 




0.20 



0.15 



0.10 



0.05 



0.00 - 




2 4 6 



10 12 14 



FIG. 2: (Color online) The density profiles of the spin-1 
^^Na condensates in the ground state for the + (solid line), 
(dashed line) and - (dash-dot lines) components respectively, 
(a) m=0, Thomas-Fermi regime; (b) m=0.2, Thomas-Fermi 
regime; (c) m=0, Tonks regime; (d) m=0.2, Tonks regime. In 
this figure the lengthes are in units of 7.4/im. 



FIG. 3: (Color online) The density profiles of the spin-1 
^^Na condensates in the ground state for the -|- (solid line), 
(dashed line) and - (dash-dot lines) components respectively. 
The interaction between atoms is enhanced such that spin in- 
teraction parameter 32 is at the same magnitude as that of 
the density interaction go- In this figure the lengthes are in 
units of 7.4 /im. 



IV. DENSITY PROFILES IN THE GROUND 
STATE 

To give a concrete example, we firstly evaluate the den- 
sity profiles of ID spinor gases in the ground state for 
*^Rb (ferromagnetic) with oq = 102aB and 02 = lOOas 
{gb is the Bohr radius) [s^- By properly choosing the 
parameters, the system may be prepared either in the 
weakly interacting regime or in the TG regime. Let us 
first consider the specific system with typical trap pa- 
rameters ojx = 0.5kHz and uj± = 50kHz for N = 2000 
atoms, in which case the effective interaction strength 
7 ~ 0.008 indicating that the system is in the weakly 
interacting regime. Fig. 1(a) and Fig. 1(b) show the 
density profiles in units of N/a for m= A4/N = and 
m= 0.2, respectively. The density profiles of -I- compo- 
nent (solid Hne) and — component (dash-dot lines) su- 
perpose on each other exactly as m= 0. When the pa- 
rameters are tuned to iVx = 10 Hz, lu± — 500kHz and 
the atomic number to iV = 50, the effective interaction 
strength 7 ~ 15 indicates that the system is in the TG 
regime. The corresponding density profiles are plotted in 
Fig. 1(c) and Fig. 1(d) for m= and m= 0.2, respec- 
tively. In the TG regime, the atoms distribute uniformly 
in more extensive area because of the strong interaction 
in the system. At the boundary the density profiles de- 
crease to zero rapidly, which reminds us the Fermi-Dirac 
statistical distribution. This imply that the density dis- 
tribution of the bosonic atoms in the TG regime behave 
like that of the Fermions. 

For the anti-ferromagnetic system we consider a con- 
densate of ^■^Na with ap = bOas and 02 = 55. las [3^ . 
In the weakly interacting regime, the trap frequencies are 



chosen as — 50Hz, uj± = lOkHz and N = 1000 so that 
7 ~ 0.001. The parameters of the system in the TG 
regime are ojx — lOHz, uj± — 2000kHz and iV = 50 with 
7 ^ 15. The density profiles are shown in the Fig. 2 in 
units of N/a for m= and m= 0.2, respectively. In the 
TG regime they exhibit similar fermionization behavior 
just as in the case of its ferromagnetic counterpart *^Rb. 
The only difference is that here the population in com- 
ponent remains completely suppressed at zero magnetic 
field according to the mean field theory. A straight expla- 
nation is that atoms prefer to be aligned anti-parallel due 
to the anti-ferromagnetic spin interaction. On the other 
hand, atoms in component would collide into pairs with 
one atom in -I- component and the other in - component 
in order to lower the spin interaction energy. This result 
is in agreement with the calculation for a spinor conden- 
sate confined in a spherically symmetric 3D harmonic 
trap, in which case the condensation occurs for + and - 
components respectively 14^]. 

For the condensates in the TG regime, the pair in- 
teraction coefficient go = N'^tt'^ /2 is clearly a constant 
irrespective of the s -wave scattering length. Should the 
interaction between atoms be enhanced by the so-called 
Feshbach resonance, the pair interaction coefficient (72 
could be increased greatly to the magnitude order of the 
coefficient go. This enable us to investigate explicitly the 
effect of the spin-spin interaction. As an example, we 
consider the system of iV = 15 ^^Na atoms and m= 0.2 
in the harmonic trap with uj^ = lOHz and uj± — 2000kHz. 
This correspond to a spin interaction parameter 172 — 500 
while go = 1110.33 if the s-wave scattering length 02 is 
enhanced to 104. Sa^. The corresponding density profiles 
are given in Fig. 3. It is shown that the density pro- 



5 



files are no longer Fermi-like and most atoms are com- 
pressed to the narrower range around the center of the 
harmonic trap. A naive explanation may be that pairs 
of atoms form singlets due to the strong spin-spin inter- 
actions and the effective interaction between the singlets 
is relatively very weak, therefore the picture of Tonks 
gas breaks down. Finally, we discuss the conditions un- 
der which the above result makes sense. For the ^^Na 
atoms discussed here, the spin exchange energy is always 
weak since C2 is two order of magnitude smaller than Cq . 
Even for the case in which 02 is enhanced to 104. 803, we 
have co/c2 = 4.14 and the theory presented in this paper 
still holds. However, for a condensate which enters the 
regime with 52 larger than go, our result is obviously not 
applicable. 



regime and the strongly interacting TG regime. The pop- 
ulation of atoms in different components depends on the 
overall magnetization and the (anti-)ferromagnetism of 
the Bose gases. When the system is in the Tonks regime, 
the density profiles show obvious Fermi-like distribution. 
However, for strong enough spin-spin interaction, we ob- 
serve apparent deviation of the density distribution from 
the Fermi-like distribution in the TG regime. 
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